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ABSTRACT 

Using an extended set of multi-epoch high resolution high signal-to-noise ratio optical 
spectra, we readdress the multiplicity properties of the 0-type stars in IC 2944 and 
in the Cen 0B2 association. We present new evidence of binarity for five objects and 
we confirm the multiple nature of another two. We derive the first orbital solutions 
for HD 100099, HD 101436 and HD 101190 and we provide additional support for 
HD 101205 being a quadruple system. The minimal spectroscopic binary fraction in our 
sample is /min = 0.57. Using numerical simulations, we show that the detection rate 
of our observational campaign is close to 90%, leaving thus little room for undetected 
spectroscopic binary systems. The statistical properties of the 0-star population in 
IC 2944 are similar, within the uncertainties, to the results obtained in the earlier 
papers in this series despite the fact that sample size effects limit the significance of 
the comparison. Using newly derived spectroscopic parallaxes, we reassess the distance 
to IC 2944 and obtained 2.3 ± 0.3 kpc, in agreement with previous studies. We also 
confirm that, as far as the O stars are concerned, the IC 2944 cluster is most likely a 
single entity. 

Key words: close - binaries: spectroscopic ~ stars: early-type - stars: individ- 
ual: (HD100099, HD101131, HD101190, HD101191, HD101205, HD101223, HD101298, 
HD101333, HD101413, HD101436, HD101545, HD308804, HD308813, CPD-62°2198) 
- open clusters and associations: individual: (Cen 0B2, IC2944) 



1 INTRODUCTION 



entities at different distances. As a result, these stars have 



Best known for its spectacular group of globules ( |Thackeray| 
1950[ [Reipurth et al.|1997| |2003j ) , IC 2944 is both an H ii re- 
gion and a rich stellar complex at the inner edge of the Ca- 
rina spiral arm. Although the stellar aggregate formed by the 
bright OB stars in the vicinity of HD 101205 (07 IIIn((f)), 



Walborn 


19731 i 


CoUinder 


(|1931|) 



1973 1 is located closer to the IC 2948 H ll region, 



( |1931| ) catalogued it as IC 2944. The same global 
group of stars was further listed as the Cen 0B2 associa- 



tion by Alter et al. ( 1970 1. The whole complex is embedded 
in the large ionization region RCW 62, that includes the 
smaller IC 2944 and IC 2948 H ll regions and extends over 
about one square-degree. Some authors (e.g., Kharchenko 
|et al.|[2005[ ) also distinguished the separate stellar cluster 
IC 2948, which is formed by the few bright stars around 
HD 101413, at ~ 12' S-E of HD 101205, and placed both 
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double entries in the stellar cluster data base WEBDA ( Mer- 
milliod||1992 1. In this paper, we follow the earlier approach 



of |Collinderf | |1931| ) and |Ardeberg fc Maurice] ( |1977| ) and we 
use the IC 2944 denomination for the complete stellar aggre- 
gate at the core of RCW 62. The homogeneity of the region 
and its distance will be discussed later in this paper. 



Thackeray & Wesselink ( 1965 1 performed the first ex- 



tensive photometric and spectroscopic study of the cluster. 
They identified 14 early B- and nine O-stars, and deduced 
a distance d ~ 2.0 ± 0.2 kpc. Several authors subsequently 
questioned the physical reality of the cluster. [Ardeberg fc| 
|Maurice| ( |1980[ ) suggested that the stellar aggregate was 
rather formed by several distinct stellar groups, with dis- 



tances ranging from 0.7 to 4 kpc. Perry & Landolt ( 1986 1 
concluded that the apparent stellar concentration of IC 2944 
resulted from the superposition of isolated early-type stars 



along the line of sight. Walborn (19871 reaffirmed, as far 



as the O-type stars are concerned, the physical reality of 
the cluster and showed that the spectroscopic parallax of 
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Table 1. Number of radial velocity (RV) measurements obtained for each object with respect to the considered spectral lines. A 
'— ' means that the corresponding line has not been measured. 
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14686 


14922 


15412 


15876 


17065 


HD 100099 


24 


24 


24 


24 


24 


24 


24 


24 


24 


14 


24 


24 


HD 101190 


20 


20 


20 


20 


20 




20 


20 


20 


10 


20 


20 


HD 101191 


12 


12 


12 


12 


12 


12 


12 


12 


12 


7 


12 


12 


HD 101223 


6 


6 


6 


6 


6 


6 


6 


6 


6 


4 


6 


6 


HD 101298 


8 


8 


4 


8 


8 


8 


8 


8 


8 


6 


8 


8 


HD 101333 


5 


5 


5 


5 


5 


5 


5 


5 


5 


3 


5 


5 


HD 101413 


10 


10 


10 


10 


10 


10 


10 


10 


10 


5 


10 


10 


HD 101436 


13 


13 


13 


13 


13 




13 


13 


13 


8 


13 


13 


HD 101545A 


7 


7 


7 


7 


7 


7 


7 


7 


7 


5 


7 


7 


HD 101545B 


7 


7 


7 


7 


7 


7 


7 


7 


7 


5 


7 


7 


HD 308804 


2 




2 


2 


2 


2 






2 




2 




HD 308813 


11 


11 


11 


11 


11 


11 


11 


11 


11 


6 


11 


11 


CPD-62°2198 


6 


6 


6 


6 


6 


5 


6 


6 


6 


4 


6 


6 



the O-type stars was in excellent agreement with the earlier 
distance proposed by Thackeray & Wesselink (19651. Po- 
larimetric observations by | Vega et al.|(|1994[) confirme d that 
some of the groups identified by Ardeberg & Maurice ( 1980 1 



showed clear differences while others could not be separated 
from each other. The more recent study of |Tovmassian et al.| 
(19981, based on the UV space telescope glazar on board 
the Mir space station (Tovmassian et al. 19881, reassessed 



the situation. Studying 185 OBA stars in an area of 8 deg 
towards Cen OB2, they confirmed that most of the O stars 
belong to IC 2944, at a distance of 2.2±0.3 kpc. They further 
identified five other OB associations with distances ranging 
from 0.85 to 6.7 kpc (see their Table 5). 

As the third paper in our series, the present work re- 
address the spectral and multiplicity properties of the O- 
type star population in Cen OB2 and IC 2944. It is orga- 
nized as follows. Sect. [2] summarizes the observing campaign 
and the data reduction. Sect. [3] discusses the individual ob- 
jects. Sect. [4] evaluates the observational biases while Sect. [5] 
describes our results and Sect. [6] summarizes our work. 



2 OBSERVATIONS AND DATA HANDLING 

This paper is based on 166 high-resolution echelle spectra 
covering most of the optical spectrum of the targets. 97 of 
those spectra were obtained from 2004 to 2007 with the 
FEROS spectrograph at the ESQ 2.2m telescope (La Silla, 
Chile) with a resolving power of 48,000. The other 69 spec- 
tra were taken from Jan to Mar 2008 using the UVES spec- 
trograph at the ESQ VLT/UT2 telescope (Paranal, Chile) 
operated in dichroic mode with the DIC2 437+760 setup 
and an 0.8" entrance slit for both arms. The resulting re- 
solving power is 50,000, thus almost identical to the feros 
one. The data reduction is similar to the one presented in 



Sana et al. ( 2009 Paper II I and the description will not be 



repeated here. The final, one-dimensional spectra were nor- 
malized to unity before merging the individual orders. The 
average signal-to-noise ratio (SNR) of the data set is well 
above 150. 

Line profile fitting was performed with a semi- 
interactive code (f itline) developed by P. Frangois at the 



Paris-Meudon Observatory and originally used in the con- 
text of chemical abundance determinations in globular clus- 
ter stars ([James et al. |2004| ) or in extremely metal-poor halo 
stars ( Frangois et al.|2007 l. This code is based on a genetic 
algorithm approach ( Charbonneau 1995 1 that mimics the 



way genetic mutations affect DNA, driving the evolution of 
species under given environmental constraints. In our case, 
lines are fitted by Gaussian profiles, each of them defined 
by four parameters: the central wavelength, the width and 
depth of the line, and the continuum value. Briefiy, the al- 
gorithm proceeds as follows: 

(i) An initial set of Gaussians is computed, giving ran- 
dom values to the four parameters. The quality of each fit 
is estimated by a calculation. 

(ii) A new "generation" of Gaussians is then calculated 
from the 20 best fits after adding random modifications to 
the initial set of parameters ("mutations"). The new set re- 
places the old one, new fits are computed and the quality of 
these fits is estimated again by using a evaluation. 

(iii) The process is iterated typically several hundred 
times until convergence to the best Gaussian fits is attained. 

(iv) A line-by-line interactive inspection of the fits is per- 
formed to remove bad fits, to correct the position of the 
continuum, or to handle line blends in SB2 systems. 

(v) The process is iterated again, and the profile param- 
eters corresponding to the best fit are adopted. 

The spectral lines considered in this process are listed in 
Table[l]together with the number of measurements available 
for each star. The equivalent width (EW) and the radial 
velocity (RV) of each line are computed based on the final fit, 
using the same rest wavelengths as adopted in other papers 
of this series. The rest wavelengths are further indicated in 
Table 2 that provides the journal of the observations and 
that lists, line-by-line, the RV measurements. 

As in other papers in the series, we adopted the spec- 



tral classification criteria of Conti & Alschuler (19711, 
|Conti| ( [T973| , |Mathys| ( pSSf and |Mathys| ( |1989l ), that 
rely on the EWs of given diagnostic lines. We used the 
usual notations: logW = logiy(A4471) - log IV(A4542), 
logW" = logiy(A4089) -logW(A4144) and logW" = 
log W^(A4388) + log W^(A4686), where the equivalent widths 
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Table 3. Diagnostic line ratios and corresponding spectral types for the studied O-type objects. Inapplicable criteria 
for specific cases are marked as n/a. We refer to text for discussion of the individual cases and to Table |8] for the 
finally adopted spectral types. Quoted error-bars give the I-ct dispersions on the measurements. 



Object 


Component 


log W 




log W" 


logVl/(A4686) 




log W" 


Sp. Type 


HD 100099 


prim 


0.305 ± 0.066 





292 ± 0.021 


n/a 


> 4 


944 ± 0.016 


09 III/I 


HD 100099 


sec 


0.660 ± 0.067 


-0 


117 ±0.054 


n/a 


> 4 


361 ± 0.062 


09.7 V 


HD 101190 


composite 


-0.369 ± 0.023 




n/a 


2.737 ±0.028 




n/a 


05.5 V 


HD 101191 


composite 


0.131 ± 0.006 





102 ± 0.032 


2.864 ± 0.005 




n/a 


08 V/III 


HD 101223 


n/a 


0.148 ±0.022 





192 ± 0.054 


2.777 ±0.026 




n/a 


08 V 


HD 101298 


n/a 


-0.200 ± 0.035 





423 ±0.143 


2.705 ± 0.037 




n/a 


06-06.5 V/HI 


HD 101333" 


n/a 


0.315 - 0.666 





118 - 0.144 


5.25 - 5.69 




n/a 


09-09.7 V/HI 


HD 101413 


prim 


0.112 ± 0.017 


-0 


038 ± 0.021 


> 2.793 ±0.013 




n/a 


08 V 


HD 101436 


prim 


-0.090 ± 0.050 





302 ± 0.230 


> 2.061 ±0.147 




n/a 


06.5-07 V 


HD 101436 


sec 


-0.170 ± 0.107 




n/a 


> 2.761 ± 0.062 




n/a 


05-06 V 


HD 101545 


A 


0.461 ± 0.041 





282 ± 0.048 


n/a 


5 


082 ±0.067 


09.5 HI 


HD 101545 


B 


0.890 ± 0.179 





234 ± 0.094 


n/a 


4 


950 ± 0.060 


09.7 III/I 


HD 308813 


prim 


0.599 ± 0.025 





015 ± 0.104 


n/a 


> 5 


378 ± 0.038 


09.5 V 


CPD -62°2198 


n/a 


0.663 ± 0.053 





205 ± 0.083 


n/a 


5 


110 ± 0.053 


09.7 III/I 



a. The nature of the adopted line profile is the dominant source of uncertainties, so that we rather quote here the 
range of values obtained by fitting either Lorentz or Gaussian profiles. 



are expressed in mA. The measurements and the correspond- 
ing spectral types are given in Table [S] and discussed in 
Sect. [3] The EWs of SB2 lines were measured at epochs of 
large separation only. Table |3] provides the actual observed 
values. For multiple systems, because of dilution by the com- 
panion continuum, the tabulated values of log M^(A4686) and 
of logW'" should be regarded as lower limits. 

In Sects. [3?T] and [3^ time series are searched for peri- 
odicities using the Fourier analysis of Heck et al. ( 1985[ ) as 



amended by Gosset et al. (20011. Quoted uncertainties on 



the period values are adopted to correspond to one tenth 
of the width of the associated peak in the periodograms. 
Whenever relevant, orbital solutions are adjusted using the 
Liege orbital solution package (e.g., Sana et al.|2006 |. 



3 O-TYPE STARS IN IC 2944 
3.1 SB2 systems 

3.1.1 HD 101131 

Close to the cluster core, HD 101131 is a known SB2 binary 



with 06.5 V((f)) and 08.5 V components. Gies et al. ( 2002 1 



published a detailed analysis of the orbital properties. They 
derived a period P of 9.6 days and a slight eccentricity e = 
0.17 ±0.03. Because HD 101131 is already well constrained, 
we did not reobserve it. 
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Figure 1. HD 101205: Hell4471, Hell AA4542, 4686 and 
He 114922 line profiles at various epochs. The right-hand column 
gives the heliocentric Julian date (HJD) at mid-exposure in for- 
mat HJD— 2 450 000. The vertical dashed lines indicate the rest 
wavelength. 



3.1.2 HD 101205 

With a combined magnitude of V = 6.45, HD 101205 is a 
visual binary with a separation of 0.36" and AV = 0.3 (|Ma- 



son et al. 1 1998 1. It is also the brightest source in the cluster 



and a known eclipsing system with a period of 2.1 d fBalona' 
1992[|Mayer et al.|1992, ptero|2007[ ). Additional photometry 
( Mayer et al. 2010 1 reveals however a more complicated hght 



curve than expected for an eclipsing binary suggesting that 
the true nature of the system is not fully understood. We 
have acquired 20 FEROS spectra from 2004 to 2006 and 15 



UVES spectra from Jan to Mar 2008. Our data reveal a clear 
multiple signature in the He I lines and an SBl motion and 
asymmetric profiles of the Hen photospheric lines (Fig. [T]). 
He nl4686 presents a mixed absorption and emission profile, 
with the main absorption component moving in phase with 
the primary and, possibly, a second fainter absorption com- 
ponent associated to the secondary. In the latter case, the 
third, overimposed emission component could be produced 
by wind-wind collision in such a short period system, a sit- 
uation similar to HD 152248 ( [Sana et al.]|2001[ ). 
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Figure 2. HD 100099: Hcll4471, Hell AA4542, 4686 and 
He 115876 line profiles at various epochs. 



The highest peak in the periodogram corresponds to a 
period of 2.8 d, a value significantly offset from the photo- 
metric period. Adopting this period, we have been able to 
fit the RVs using a circular orbital solution (as expected for 
such a short period system), yet the residuals are widely 
spread: 18 kms~^. This is significantly larger than the ex- 
pected RV accuracy (< 5 kms~^). This suggests either a 
systematic error in the measurements (e.g., the presence of a 
third spectroscopic component), or an additional variability 
signal. Naturally, the presence of the speckle companion im- 
plies that a third Hght spectral component must be present, 
and line blending would tend to reduce the measured veloc- 
ity excursions, ft is unclear whether it would however bias 
the periodicity of the signal. Finally, a multi-periodic analy- 
sis of the time series reveals however no extra-periodic signal 
on the time scales covered by our campaign, beside the 2.8 d 
period. 

Our data set is unfortunately insufficient to solve the 
discrepancy between the photometric and spectroscopic pe- 



riods. As noted by Mayer et al. ( 2010 1, the two visual compo- 



nents of HD 101205 might both be multiple systems. A ded- 
icated monitoring is needed to elucidate the nature of this 
object, which could be a new Trapezium system. The mul- 
tiple nature of HD 101205 is, however, no longer in doubt. 

3.1.3 HD 100099 

At about 1° SW of HD 101205, H D 100099 is associated to 
IC 2944 in the earlier works (e.g. , [Thackeray fc Wesselink| 
1965 1 although Tovmassian et al. ( 1998 1 preferred to assign 



it to the Cen-Cru OB2.7 association. Based on the large 
angular separation of HD 100099 from the cluster center. 



Baumgardt et al. (20001 assigned to it a low membership 
probability but noted that its proper motion was still com- 
patible with the kinematic properties of IC 2944. 

HD 100099 was classified 09 III by [Garrison eFal] 



Table 4. HD 100099: best-fit orbital solution based on the 
He 115876 RV measurements. T (in HJD-2 450 000) is the time 
of periastron passage and is adopted as if) = 0.0 in Fig. |4] 
Quoted uncertainties correspond to 1-it error-bars. 



4521,779 




He 115876 


4514,883 


p(d) 


21.5585 ±0.0330 


4493,830 




0.517 ±0.019 






305.5 ±0.9 


3864.501 


T 


2995.654 ± 0.078 


3863.489 


7i (kms"-*-) 


16.5 ± 2.1 




72 (kms~^) 


-6.0 ± 2.6 


3862.485 


Kx (kms-i) 


138.4 ± 5.7 


3861.510 


K2 (kms-i) 


174.8 ± 7.3 




M2/M1 


0.792 ±0.013 


3860.551 


ai sini (R©) 


50.46 ± 2.20 


3511-507 


02 sini (R0) 


63.72 ± 2.78 




Ml sin^ i (Mq) 


24.0 ± 3.1 


3135.461 


M2 sin^ i (Mq) 


19.0 ± 2.4 




rms (kms~^) 


2.4 
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Figure 3. HD 100099: periodogram (upper panel) and spectral 
window (lower panel) computed from the Heli5876 primary RV 
measurements. 




Figure 4. HD 100099: primary (filled symbols) and secondary 
(open symbols) RVs together with the best-fit RV curves. The 
crosses indicate the data not taken into account in the fit because 
components could not be securely separated. 
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Figure 5. HD 101436: Hcll4471, Hell AA4542, 4686 and 
He 115876 line profiles at various epochs. 



( 1977[ ) and by |Herbst||l975| ). We collected 24 spectra of this 
star that clearly revealed an SB2 system with two late O- 
type components (Fig.[2|. Clear RV variations are seen from 
one night to the other, suggesting a relatively short period. 
Both Hell and Siiv lines are definitely present in the sec- 
ondary spectrum. Disentangling the secondary signature in 
the Heiil4542 region is more difficult because of the pres- 
ence of N III AA4534-4545 in the primary spectrum that has 
a similar strength as the secondary Helll4542 line. Reliable 
separation is achieved on a couple of spectra only. Our best 
estimate for the observed secondary He Ill4542 EW is O.OSA, 
and we used this value to derive the spectral type. Based 
on our measurements, the spectral type estimates are 09 
and 09.7 for the primary and secondary stars respectively, 
with the 08.5 and 09.5 spectral sub-types well within the 
respective uncertainties. The Siivl4089 over Hell4144 line 
ratio, W" , indicates that the primary is likely a giant, with 
the supergiant class within the uncertainties, and that the 
secondary is likely a main sequence star. Given the object's 
magnitude {V = 8.09), a supergiant primary would place 
the object at 4 to 5 kpc and a class HI would place it at 
about 3 kpc. While HD 100099 is thus potentially a mem- 
ber of a different association (Tovmassian et al. 19981, its 



systemic velocity is compatible with the RVs of the other O 
stars in IC 2944. 

A detailed study of the orbital solution is beyond the 
scope of this paper. Here we limit ourselves to a preliminary 
analysis based on the Hell5876 line. Our Fourier analysis 
indicates a most probable period close to 21.57 d, although 
we note that the corresponding peak in the periodogram is 
favoured by the sampling (Fig. [sl . Best-fit parameters and 
RV curves are displayed in Tablejljand in Fig. [4] respectively. 
While the fit is of good quality, our data do not sample the 
maximum separation epoch, leaving some uncertainties both 
on the RV curve semi-amplitudes Ki and K2 and on the 
eccentricity. 



Figure 6. HD 101436: primary (filled symbols) and secondary 
(open symbols) RVs together with the best-fit RV curves of Ta- 
ble H 



3.14 HD 101436 



HD 101436 (06.5 V, |Walborn|19"73| was flagged as a p roba- 
ble RV variable star by [Thackeray &: Wesseliiik| ( |1965| , with 
significant RV changes on a time scale of a week. Two ad- 



ditional RV measurements by Humphreys ( 1973 1 and Conti 
(19771 confirmed a ARV amplitude of Ri40 kms"\ 



et al. 



We obtained 13 additional spectra covering time scales from 
days to years that reveal, for the first time, a shallow sec- 
ondary signature in the He I, Hen and Silv lines (Fig. [5|. 
Fitted with a single Gaussian, the narrow (primary) com- 
ponent presents a clear and smooth increase over our 5-day 
campaign in May 2006, rising by about 40 kms~^ from —20 
to +20 kms~^ approximately, although slightly depending 
on the line considered. Yet, RV measurements are hampered 
by a degeneracy in the determination of the line profiles. 
Two different hypotheses allow us to fit the data. Either the 
secondary contributes to the blended profile as a faint and 
narrow component or it contributes as a broad and shallow 
one. In both cases however, the values of the Conti's cri- 
teria are left unaffected within the uncertainties, with the 
primary being an 07 star and the secondary, an 06.5 star 
with half sub-type uncertainties on both measurements. Be- 
cause the Hell4144 line is very weak, the W" criterion is 
unsuited and we relied on the W^(A4686) measurements to 
estimate the luminosity class. Depending on the adopted 
width for the secondary, the He nl4686 line is dominated ei- 
ther by the primary or by the secondary star. The dominat- 
ing component displays a strong absorption, typical of dwarf 
stars. The other component shows then a measured value of 
log W(A4686) ^ 2.2 - 2.3. Because the hne is diluted by 
the companion continuum, this measurement provides only 
a lower limit on the true strength of the Helll4686. It would 
correspond to a main sequence star if the light ratio is 1:3 or 
more. Given the system magnitude, both stars are further 
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Table 5. HD 101436: best-fit orbital solution using RVs 
based on Gaussian profiles (Col. 2) and on atmosphere model 
profiles (Col. 3). T (in HJD-2 450 000) is the time of peri- 
astron passage. Quoted uncertainties correspond to Icr error- 
bars. 



Parameter 


Gaussian profiles 


Atmosphere model profiles 


P(d) 


37.37 ±0.14 


37.37 ±0.14 


e 


0.11 ± 0.01 


0.11 ±0.02 


oj n 


274 ± 12 


296 ± 11 


T 


3190.77 ± 1.14 


3193.04 ± 1.01 


M1/M2 


0.52 ±0.02 


0.79 ± 0.03 


71 (kms^^) 


6.0 ±0.8 


3.2 ±0.8 


72 (kms^-') 


-8.7 ± 0.5 


-2.2 ± 0.6 


Ki (kms-i) 


46.9 ± 1.4 


43.4 ± 1.1 


K2 (kms-i) 


24.3 ±0.7 


34.4 ± 0.9 


rms (kms~^) 


1.3 


1.4 



expected to be main sequence stars if located in IC 2944 and 
we adopt 07 V+06.5 V as a tentative spectral classification. 

To double check the results of the RV genetic fitting 
algorithm, we also developed a global minimization ap- 
proach. We adjusted the profiles of the He: AA4026, 4471, 
5876 and of the Hen AA4200, 4542, 4686 lines simultane- 
ously on all the 13 spectra using the following assumptions: 

(i) the profile of each line is formed by a blend of two Gaus- 
sian components representing the contribution of each star, 

(ii) the RV of each star is identical for all the lines at a given 
epoch, (iii) the amplitudes and full widths at half maximum 
(FWIfMs) of the Gaussians are identical for a given line and 
a given stellar component at all epochs. As a result, this 
approach definitely favours the second option of a narrow 
primary and a broad and shallow secondary. 

As a final step, we replaced the Gaussian profiles by 
more realistic profiles. The latter were computed using 



FASTWIND (Puis et al. 20051 for a range of parameters 



typical of mid-O stars ( Martins et al.|2005 l, then convolved 



with a rotational profile. Given the high spectral resolution 
of our data, we neglected the effect of the instrumental pro- 
file. Dilution by the companion continuum was included as a 
free parameter. To increase the realism of the approach, we 
further required that each line of a given component uses 
the same dilution factor and the same projected rotation 
rate. The best fit was obtained by using atmosphere mod- 
els corresponding to the physical parameters of an 06 V 
primary with vaini « 70 kms~^ and an 06.5 V secondary 
with vsini « 200 kms~^. 

While the overall line profiles are well reproduced in 
both approaches, the obtained RVs show systematic differ- 
ences. On average, the RVs associated to the narrow com- 
ponent are 3 kms~^ smaller when using the FASTWIND 
profiles and the RVs of the broad component are 8 kms^^ 
larger. This has a significant impact on the respective am- 
plitudes of the RV signal as well as on their ratio (Table 

Despite these uncertainties, the global x minimization 
improves the disentangling performances compared to the 
genetic algorithm described in Sect. [2j allowing us to con- 
strain some of the orbital parameters. For both RV sets, a 
Fourier analysis indicates two most probable periods of 12.9 
and 37.3 d, the second one presenting the highest peak in 
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Figure 7. HD 101190: Hell4471, Hell AA4542, 4686 and 
He 115876 line profiles at various epochs. 



the periodograms. A preliminary orbital solution analysis 
also favours the longer period. In both cases, slightly eccen- 
tric orbits provide a better fit compared to circular orbits. 
All in all, the orbital solutions built on both RV sets yield 
parameters in good agreement, except for the amplitude of 
the secondary RV curve and for the mass-ratio. In particu- 
lar, the RV measurements based on the more realistic FAST- 
WIND profiles yield a mass-ratio closer to unity, which is in 
better agreement with our spectral type estimates. 

Fig. [6] displays the best-fit RV curves and Table [5] lists 
the corresponding parameters. We note that we do not con- 
strain the maximum RV separation, which is a source of 
uncertainty for the eccentricity and for the semi-amplitude 
of the RV curves. While we consider that the period is well 
constrained and that the system has likely a limited or no 
eccentricity, we regard other parameters as preliminary. 

Under the adopted hypotheses, the component that we 
have referred to as the secondary appears to be the hotter, 
most massive component. Given the low amplitude of the 
RV curve, the minimum masses are very small, with values 
slightly below 1 Mq, suggesting an inclination of the order 
of 20° . This implies a deprojected rotational velocity for the 
fast rotating component of about 600 kms^^. One possi- 
ble channel to produce such a fast rotator is a recent mass 
transfer from the primary to the secondary that would have 
increased the secondary mass and spin it up to the critical 
rotation rate. If this scenario is correct, HD 101436 could 
turn out to be a very interesting system to constrain evo- 
lutionary scenarios of massive binaries. More data and, in 
particular a better coverage of the extrema of the RV curves, 
and the use of a disentangling algorithm to recover the line 
profiles independently of any assumption on their shapes are 
the needed ingredients to verify our preliminary analysis. 
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Figure 8. HD 101190: periodogram (upper panel) and spec- 
tral window (lower panel) computed using the average of the 
He I AA4026, 4471, 4922, 5876 line RVs. 
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Figure 9. HD 101190: best-fit RV curve for the He I (upper 
panel) and He II (lower panel) data sets. 



Table 6. HD 101190: best-fit orbital solutions based on the 
He I and Hen data sets. T (in HJD-2 450 000) is the time of 
periastron passage and is adopted as </> = 0.0 in Fig.|9] Quoted 
uncertainties correspond to Icr error-bars. 



Nv AA4604-4620 absorptions and of the Siiv AA4089-4116 
emissions is not compatible with the 06 V((f)) type reported 



Parameter 


He I lines 


He II lines 


P(d) 


6.0466 ± 0.0026 


6.0473 ± 0.0026 


e 


0.29 ±0.04 


0.36 ±0.07 




117.8 ±4.8 


127.5 ±6.2 


T 


2999.17 ± 0.09 


2999.13 ±0.11 


72 (kms^^) 


-3.8 ±0.7 


8.2 ±0.7 


K2 (kms-1) 


34.4 ± 1.4 


20.8 ± 1.2 


rms (kms^^) 


2.5 


2.1 



3.1.5 HD 101190 

To the North of the cluster, HD 101190 (06 V((f)), 



born 1973 \ is reported as RV variable star ( Ardeberg & 



Wal- 



Maurice 19771. Photometric monitoring by Balona (19921 
also revealed significant variability with a rms amplitude of 
7 millimag although no period has been identified so far. 
We collected 20 spectra over a 4 year timebase that shows 
relatively narrow lines, consistent with the rotational veloc- 



ities deduced by Penny ( 1996 \ and Howarth et al. ( 1997 1 
(?;sini=89 and 88 kms~^respectively). Our data further re- 
veal clear variability of all the considered lines with asym- 
metric profiles seen at several epochs (Fig.[7|. The He I lines 
show a larger RV amplitude than the He ll lines but no sec- 
ondary signature could be reliably disentangled by Gaussian 
fitting. Our classification criteria assign an 05.5 V((f)) spec- 
tral type to the composite spectra. 

Inspection of the spectra reveals several interesting fea- 
tures. All the spectral lines in the HD 101190 spectrum dis- 
play a clearly correlated motion, except the Nlll AA4634- 
4640 emission lines and the Nv AA4604-4620 absorption 
lines that show constant velocity (o"rv4604,4620 ~ 3 kms^^). 
On two occasions, the Silv AA4089, 4116 lines are seen 
separated into a blue-shifted absorption and a red-shifted 
emission, the latter emissions displaying RVs consistent 
with those of the Nlll and Nv lines. The presence of the 



by Walborn ( 1973 1, nor is it compatible with the average of 
the Hell4471 over Helll4542 ratio measured from our data. 
This suggests that the HD 101190 spectrum has two compo- 
nents: a very early O-type star that shows no RV variations 
and that is responsible of the N v absorptions and Si iv emis- 
sions, and a later O-type star which carries the RV signal. 

Because we cannot detect the Nivl4058 line and be- 
cause the Nv and Silv lines are relatively faint, we assign 
the hotter star an 04 V((f+)) type, although we considered 
the 03.5 sub-type within the uncertainties. Given the ex- 
pected visible flux ratio between an 04 V star and a less 
massive 07 component (« 2.75, [Martins et al.||2005 ), one 
finds that an 07-07.5 V star, with « 1, can reproduce 
the observed W' ratio of the composite spectrum. While our 
data do not allow us to prove that the 04 and 07 compo- 
nents are physically bound, we note that the chance align- 
ment of two such stars is very small. 

In an attempt to constrain the orbital properties of the 
system, we performed a Fourier analysis of the RV series as- 
sociated with the individual He I and He II lines. All the data 
sets consistently yield an orbital period of 6.047 d, which 
is indeed the dominant alias in the periodogram (Fig. [8|. 
We computed SBl orbital solutions both using the RVs 
from the individual lines and by combining separately the 
He I AA4026, 4471, 4922, 5876 and the Hell AA4200, 4542, 
4686 RV measurements. Again we obtained very similar re- 
sults. Table |6] and Fig. |9] give the best fit orbital solution 
for the He I and Hell data sets and show the corresponding 
RV curves. While we consider that the orbital period is re- 
liable, we note that the semi-amplitude of the RV curve K, 
is likely a lower limit as it is affected by the blend with the 
04 signature. 

Finally, assuming that the secondary is not a very late 
O-star, an hypothesis supported by the RV signatures in 
the He ll lines, the binary mass-ratio AI2 /M\ is likely larger 
than 0.3. This would imply a lower limit on the primary RV 
curve semi-amplitude: Ki > 10 kms~^, so that one would 
have expected to see some variations of the primary lines. 
While this argues against a direct bound between the 04 and 
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Figure 10. HD 101413: Hell4471, Hon AA4542, 4686 and 
He 115876 line profiles at various epochs. 



Figure 11. HD 101191: Hell4471, Hell AA4542, 4686 and 
He 115876 line profiles at various epochs. 



the 07 star, we note that a mass ratio a factor of two lower 
than assumed here would solve this apparent discrepancy. 
Alternatively, the system may be triple, with a single hot 
star and an SBl cooler star. More data are definitely needed 
to solve this issue. 



3.1.6 HD 101413 



HD 101413 (08 V, Walborn 19731 is the visual compan 



ion of HD 101436, with a separation of 27.8" (Mason et al 
|1998[ ). [Thackeray fc Wesselink] ( |1965| ) first reported large 
RV variations {ARV ~ 100 kms~ ) based on 3 spectra col- 
lected over a year. lUE and Fuse rot ational velocities were 
measured in the range 90-102 kms"^jPennyll996l Howarth 



et al.|1997| [Penny fc Gics 2009) and probably relate to the 



primary or to the composite spectrum. We obtained three 
FEROS spectra in May 2005, two in May 2006 and a set of 
5 UVES spectra spread over 3 months early 2008. The 2005 
and some of the 2008 observations revealed the signature of 
the secondary component for the first time, although only in 



the He I lines (Fig. 10 1. A limited RV difference is seen on a 



day-to-day basis, but the primary and secondary smoothly 
varied from Jan to Mar 2008, suggesting a period of 3 to 6 
months. We note that we never observed such an extreme 
primary RV as the one reported by [Thackeray fc Wesselink[ 
(1965[) in Apr 1957 {RVprim = -81 kms"'). Using the few 



well separated spectra, we confirm the primary 08 V clas- 
sification. The companion spectrum shows neither He 11 nor 
Si IV lines and we assign it a B2-3 V classification. While we 
are missing a complete orbital coverage, empirical determi- 
nation of the mass-ratio can still be obtained. Based on the 
Hel lines, we obtained M2/M1 — 0.17 ± 0.01, which rather 
suggests a mid- to late-B type for the secondary. 



3.2 SBl systems 
3.2.1 HD 101191 



HD 101191 (08 V((n)), |Walborn|l973[ ) is located at 1.4' SW 
of HD 101205, in the core of the cluster and has been scarcely 
studied so far. We obtained 12 high-resolution spectra that 



reveal a clear SBl, long period signature (Fig. 111. Our dif- 



ferent epochs indeed reveal RV variations with a peak-to- 
peak amplitude of ~30 kms~^ for the Hell lines and of 
~20 kms^^ for the Hel lines. No significant night to night 
variability is observed, nor within the 3-month coverage of 
our UVES data. While we cannot rule out an alias at 20 d, 
such a period value is unlikely because the 3- and 5-day 
FEROS campaigns in May 2005 and May 2006 reveal no RV 
variation at all. This rather suggests a long orbital period, of 
several times the 3-month timebase of our UVES campaign. 
From our data, we can also place a relatively firm upper 
limit on the orbital period of about 4 yr. As the measured 
RVs display a relatively limited amplitude, we might have 
missed the epoch of maximum separation, keeping open the 
possibility that one can still observe the signatures of the 
two components. We derived an 08 V/III type for the com- 
posite spectrum. Because the He nl4686 line is strongly seen 
in absorption, we adopt the 08 V spectral sub-type as our 
final classification of the HD 101191 composite spectrum. 

3.2.2 HD 308813 

At 4.2' NW of HD 101205, HD 308813 was classified 09.5 V 



by[Schild|([l970|) and BO V byJArdeberg fc Maurice[ ( [1977 1. 
Thackeray & Wesselink ('1965) observed it three times over 
a 2 yr time-scale and reported IS.RV ~ 28 kms~^. Based 
on discrepant Hell RVs compared to the Hel ones, they 
proposed that HD 308813 is an SB2 candidate. [Huang fc| 
[Giesj ( [2006a[ ) also collected three observations over 4 days 
and reported RVs from 21.9 to —17.6 kms~^, confirm- 
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Figure 12. HD 308813: Hell4471, Hon AA4542, 4686 and 
He 115876 line profiles at various epochs. 



ing its SBl nature. They further measured a projected 
rotational velocity v sin i of 196 kms-^ We collected 11 
spectra of HD 308813, covering time scales from days to 
years. HD 308813 presents broad, somewhat variable lines 
(Fig. |12[ ). We measured RVs using both Gaussian profile 
and purely rotationally broadened profiles and obtained very 
comparable results. Clear RV shifts (peak-to-peak difference 
over 40 kms"^) are seen, although these are not happening 
as quickly as reported by Huang & Gies (2006a I. All the 



lines follow the same orbital motion, so that we cannot con- 
firm the SB2 nature of the object. We however note that 
some lines display systematic shifts in their RVs compared 
to the average. Given our data set, we cannot distinguish 
between a very-short period system (P < 2 d) and a longer 
one (P > 20 d). With a clear Helll4542 line, HD 308813 is 
not a B star but a late O star. We thus confirm the 09.5 V 
classification of |Schild| ( |l970[ ). 
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Figure 13. HD 101298: Hell4471, Hell AA4542, 4686 and 
He 115876 line profiles at various epochs. 
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Figure 14. HD 101223: He 114471, Hell AA4542, 4686 and 
He 115876 line profiles at various epochs. 



3.3 Presumably single stars 

3.3.1 HD 101298 



The object was observed spectroscopically b y [Thackeray 
fc Wesselink| (|1965[), [Conti et aX] (|1977[) and [Ardeberg 



Maurice (19771. The latter reported possible RV variations 



within their four measurements. A 6-night photometric mon- 



itoring ( Balona 1992 1 however showed a constant luminosity. 



Similarly, no significant RV difference is observed in between 
the eight epochs of our campaign, which again suggests that 



the star is probably single (Fig. 13 1. Spectral type estimates 
point toward an 06-6.5 star. The W" criterion is unsuited 
for such an early star as the Hell4144 line has become very 
weak. With log W(A4686) « 2.7, the Heiil4686 line suggests 
a giant classification. Being slightly variable while the other 
lines remain constant, the Helll4686 line further suggests 
the presence of a significant wind, another hint of a slightly 



more evolved object. Following Walborn (19731, we adopt 



an 06 lll((f)) classification. Note however that a class V 
would be in better agreement with the overall distance to 



the cluster (see Sect. 5.1 1 



3.3.2 HD 101223 



Ardeberg & Maurice (19771 reported discrepant velocity 
measurements (ranging from —29 to -1-4 kms~^) and asym- 
metrical line profile, suspecting HD 101223 to be a spec- 
troscopic binary. Observed at six epochs over five years 
(Fig. |14[ ), the HD 101223 spectrum remains constant with a 
Icr RV dispersion of the order of 1 kms~^. The classification 
criteria indicate an 08 V/llI star. Because the Helll4686 
line is strongly seen in absorption, we prefer the 08 V clas- 
sification and we consider the star as single. 
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Figure 15. CPD-62°2198: Hell4471, Hcll AA4542, 4686 and 
He 115876 lino profiles at various epochs. 



Table 7. Binary detection probability for the time sampling asso- 
ciated with different objects (Col. 1) and for various period ranges 
(Cols. 2 to 5). 



Time sampling 


Short 


Intcrmed. 


Long 


All 




[2-lOd] 


[10-365d] 


[365-3000d] 


[2-3000d] 


HD 100099 


0.995 


0.920 


0.631 


0.902 


HD 101131 


0.995 


0.970 


0.813 


0.954 


HD 101190 


0.996 


0.944 


0.730 


0.930 


HD 101191 


0.992 


0.906 


0.614 


0.894 


HD 101205 


0.995 


0.893 


0.560 


0.884 


HD 101223 


0.988 


0.864 


0.607 


0.877 


HD 101298 


0.994 


0.905 


0.688 


0.910 


HD 101333 


0.985 


0.759 


0.324 


0.789 


HD 101413 


0.992 


0.900 


0.530 


0.876 


HD 101436 


0.995 


0.915 


0.568 


0.891 


HD 101545A 


0.991 


0.783 


0.376 


0.810 


HD 101545B 


0.992 


0.785 


0.341 


0.805 


HD 308813 


0.992 


0.886 


0.463 


0.859 


CPD-62°2198 


0.987 


0.839 


0.511 


0.851 



3.3.3 CPD~62°2U 



3.3.5 HD 101545A & B 



At about 30' NW of the cluster center, CPD-62°2198 
presents a significantly different RV compared to other sin- 
gle stars (from about —6 to —12 kms~^ depending on the 
considered lines). Yet, we could not detect any significant 
RV variations within our 6 observations (Fig. |15[) no r by 
comparison with older data by Reed & Kuhna ( 1997 1 and 



Huang & Gies (2006a I . Classified as 09.5 lb by Ardeberg 
fc Maurice] ( |1977[ ), our Hell4471 over Helll4542 ratio in- 
dicates an 09.7 star, with the 09.5 sub-type at la. The 
iy(A4089)/lV(A4144) ratio indicates a class III star while 
the log W'" criterion rather points to a class I star with 



class III within la. Huang & Gies (2006b I obtained T^s — 



29132 ±375 K and logp = 3.556 ±0.039 for CPD-62°2198, 
in good agreement with values expected for an 09.7 III star. 
We consider CPD— 62°2198 to be likely single and we finally 
adopt the 09.7 HI classification. 



3.3.i HD 101333 

Classified BO HI 
39 kms"^ 



HD 101333 is a slow rotator {vsmi 



Huang & Gies 2006a) in the Southern part of the 



cluster. Thackeray fc Wesselink| ( 1965 I and Huang & Gies 
( |2006a| ) both reported a set of three RV measurements that 
show no internal discrepancy but these three sets differ by 
15 kms~^ from each other. We obtained three feros spectra 
in May 2006 and another two UVES spectra in January and 
March 2008. Several lines (e.g. He: AA4144, 4388) display 
a clear Lorentz profile while others have a well pronounced 
Gaussian shape. The UVES data are about 5 kms~^ larger 
than the feros measurements. Given the slow rotation, this 
is definitely significant although we cannot decide between 
a low-amplitude photospheric effect and a large mass-ratio 
and/or low inclination binary scenario. The obtained spec- 
tral classification is somewhat dependent on the assumed 
line profile, but points to a late-type O star. Comparing 
HD 101333 with the spectral standards t Sco (BO. 2 V) and 



HD 98028 (09 V) from |Walborn fc Fitzpatrickl ( [1990| ) and 
to BD-13°4930 (09.5 V) from |Paper II| (Fig.|16]), we finally 
adopt an 09.5 V classification. 



HD 101545 is a visual binary the components of which are 
separated by a few arcsecs, the component A being to the 



NE of component B. Mason et al. ( 1998 1 reported a sepa- 



ration of 2.57 and respective V mag of 6.9 and 7.4. Both 



components are late O-type stars (Fig. 17 1. None of them 



display any significant RV variation in their respective time 
series (7 spectra over 5 yr for each object of the pair). The 
best spectral type estimate yields 09.5 III/I and 09.7III/I 
for the A and B components. Given the stars' brightness, we 
finally consider both stars to be giants and not supergiants. 

3.4 Other stars in IC 2944 

3.4.1 HD 



Reported as 09.5 V by Ardeberg fc Maurice (19771, we 
observed HD 308804 twice with a 2 d separation. The 
HD 308804 spectrum contains neither Hell nor Siiv lines 
and it is thus not the spectrum of an O star. While [Arde-| 
berg & Maurice ( 1977 1 mentioned double or asymmetric 
lines, our two spectra are single-lined. Yet, they display a 
RV difference of ~ 10 kms~^, although the first spectrum 
has a much lower SNR compared to the second one. |Huang| 
fc Gies (2006a I also observed this star three times over a 



four-night period and reported no RV change. Their mea- 
sured RV is however 10 and 20 kms~^ more negative than 
our two measurements, providing further support to the bi- 
nary hypothesis. By comparison with the atlas of |Walborn fc| 
Fitzpatrick ( 19901, we adopt a B3 V type for the composite 



spectrum of HD 308804. 



4 OBSERVATIONAL BIASES 

To estimate the impact of the observational biases on the de- 
tection probability of a binary with a given primary mass, 
we follow the approach described in [Paper II] Briefly, we 
first draw a series of A'^ SB2 systems adopting the following 
rules. The primary mass is fixed to the theoretical value of 



Martins et al. ( 2005 1 according to the primary spectral type 
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Figure 16. HD 101333: 4000-4700 A spectral range compared to the spectroscopic standards r Sco and HD 98028 and to BD-13°4930. 
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Figure 17. HD 101545A (left) and HD 101545B (right): Hcll4471, Hen AA4542, 4686 and Hell5876 line profiles at various epochs. 



of the considered object. The probabihty density function 
(PDF) of logP (d) is taken to be uniform in the interval 
[0.3,1.0] and [1.0,3.5]. 50% of the systems are spread in the 
first interval while the rest is distributed over the second in- 
terval. The eccentricity and mass-ratio PDFs are taken to be 
uniform in the range [0.0,0.9] and [0.1,1.0] respectively, with 
the additional constraint that the separation at periastron 
is larger than 20 R0 to avoid contact. Finally, the orbital 
configuration of each system is drawn assuming a random 
orientation in space and a random time of periastron pas- 
sage. 

We then applied the time sampling corresponding to the 
observations of a given object to the RV curve of each of the 
A*' SB2 binaries and we computed what would have been 
the largest RV difference (Ai?y) recorded. The detection 
probability Pdotcct is then given by the ratio of the systems 
that display ARV larger than a given threshold C: 



NjARV > C) 
N ■ 



(1) 



In the results presented in Table [7[ we adopted C — 
20 kms~^ and we computed Pdotoct for various ranges of the 
orbital periods. Our overall detection probability is better 
than 0.8 for all the objects, with the main uncertainties com- 
ing from the long period range. Based on those results one 
can exclude that the objects reported as single (Sect. 3.3 1 



are short period binaries (P < 10 d). Some of them could 
still be undetected long-period binaries. Yet, given that such 
systems are expected to be rare ( [Sana fc Evansl|2010[ ), the 
chance of having at least one undetected system among the 
six presumably single objects is 0.56. 

Fig. [18] displays the probability density function of de- 
tecting a given number of binaries assuming all the objects 
in IC 2944 are multiple systems. It reveals that the probabil- 
ity to detect 8 binaries or less is 0.002, allowing thus to reject 
the null hypothesis of /true = 1.0. This plot also indicates 
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Figure 18. Distribution of the number of detected SB systems 
as obtained from 10 000 realizations of a cluster with 14 O-typc 
binaries. 

that the most hkely binary detection rate of our campaign 
is close to 0.9. 

As discussed in |Paper II[ our simulations neglect the ef- 
fect of line blending, that decreases the detection probability 
of near equal flux systems and/or of systems with large pro- 
jected rotational velocities. Among our sample, HD 308813 
displays the broadest lines («sinj ~ 200kms~^), yet RV 
variations with a peak-to-peak amplitude of ~ 40 km s~^ are 
clearly detected. All the presumably single stars discussed 
in Sect. [33] have significantly narrower line profiles so that 
even small RV variations in an equal-fiux system would in- 
duce significant changes in these line profiles. Finally, we 
note that the main limitations of our simulations lay in the 
adopted distributions of orbital parameters. While our as- 
sumptions are in line with earlier results from our campaign 
as well as with the preliminary results of [Sana fc Evans| 
( |2010[ ), it is unknown whether the multiplicity properties of 
O stars are homogeneous or show sizeable variations from 
one cluster to another. 



Figure 19. Dispersion of the Hell5876 RV measurements versus 
their mean RV. Circles around the a — fi points indicate the bi- 
naries. The average systemic velocities have been used instead of 
the mean RV for the binaries with orbital solutions : HD 100099, 
HD 101131, HD 101190 and HD 101436. 



Table 9. Color excess and spectroscopic distance modulus 
(DM) computed from the V and K band magnitudes. 



Object 


E{B - V) 




DMk 


HD 100099 


0.36 


12.57 ± 0.26 


12.31 ± 0.07 


HD 101131 


0.27 


11.68 ±0.40 


11.40 ±0.11 


HD 101190 


0.37 


12.08 ±0.50 


12.12 ±0.21 


HD 101205 


0.30 


11.43 ±0.73 


11.08 ±0.10 


HD 101413 


0.30 


11.86 ±0.28 


11.58 ±0.13 


HD 101436 


0.32 


12.14 ± 0.50 


12.00 ±0.11 


HD 101191 


0.28 


12.03 ±0.14 


11.79 ±0.14 


HD 308804 


0.40 


11.73 ±0.55 


10.82 ±0.05 


HD 308813 


0.29 


12.28 ±0.16 


12.17 ± 0.09 


HD 101223 


0.43 


11.77± 0.15 


11.63 ±0.14 


HD 101298 


0.31 


12.86 ±0.09 


12.57 ± 0.08 


HD 101333 


0.31 


12.02 ±0.16 


11.75 ±0.09 


HD 101545AB 


0.23 


12.16 ±0.32 


11.40 ±0.05 


CPD-62°2198 


0.39 


14.46 ±0.07 


13.05 ±0.05 



5 DISCUSSION 

5.1 Distance and membership 

As already discussed in the Introduction, several authors 
questioned the single nature of the IC 2944 complex and 
preferred to rather consider several OB associations along 
the line of sight. Using our refined spectral classifications, 
we recompute the distance modulus (DM) of each target 
in our sample using both the V and K band magnitudes 
(Table [9]|. In this exercise, we adopted the theoretical ab- 
solute magnitudes, colors and bolometric corrections from 



Fitzpatrick ( 1999 1 



Martins & Plez ( 2006 1 and the interstellar extinction law of 



The revision of the spectral classification and more 



particularly of the multiplicity status of the components 
leads to significant changes in the spectroscopic parallax 
of the objects. The most critical examples are probably 
the cases of HD 101191 and HD 101545 which are now lo- 
cated at a distance compatible with the other objects in 
IC 2944. The largest deviations are shown by HD 308804 
and CPD — 62°2198 and these stars might respectively be 
foreground and background objects. HD 101298 also seems 
slightly more distant. As noted earlier, adopting a main- 
sequence classification instead of a giant one would bring 
the star in perfect agreement with the cluster mean dis- 
tance. Ignoring these three stars as well as HD 101205 due to 
its uncertain number of components, the estimated distance 



14 H. Sana et al. 



Table 8. Final spectroscopic classification and multiplicity properties of the studied O stars in IC 2944. An 'O' 
appended to the SBl/2 flag means that an orbital solution is available; an 'E', that the system displays eclipses in 
addition to its SB signature. 



Object 


Mult. 




Sp. Type 




P 


e 






Previous works 


This work 




(d) 




SB2 systems 


HD 100099 


SB20 


09 III 


09 ni+09.7 V 


0.792 


21.56 


0.52 


HD 101131 


SB20 


06.5 V((f))+08.5 V n/a 


0.56 


9.65 


0.16 


HD 101190 


SB20 


05.5-06 V({f)) 


04 V{(f))+07 V 


0.3-0.5 


6.05 


ft! 0.3 


HD 101205 


SB20E 


07 IIIn((f)) 


n/a 


undef. 


2.1-2.8: 


0.0: 


HD 101413 


SB2 


08 V 


08 V+B3: V 


^ 0.2 


1 - 2 X 10^ 


undef. 


HD 101436 


SB2 


06.5 V 


06.5 V+07 V 


0.52-0.79" 


37.37 


0.11 


SBl systems 


HD 101191 


SBl 


08 V{(n)) 


08 V 


undef. 


10^ - 10^ 


undef. 


HD 308804 


SBl 


09.5 V 


B3 V 


undef. 


> 5 


undef. 


HD 308813 


SBl 


09.5 V-BO V 


09.5 V 


undef. 


undef. 


undef. 


Single stars 


HD 101223 


sgl. 


08 V{(f)) 


08 V 


n/a 


n/a 


n/a 


HD 101298 


sgl. 


06 ni((f)) 


06 ni((f)) 


n/a 


n/a 


n/a 


HD 101333 


sgl. 


O9.5-B0 III 


09.5 V 


n/a 


n/a 


n/a 


HD 101545A 


sgl. 


O9.5-B0.5 lb 


09.5 HI 


n/a 


n/a 


n/a 


HD 101545B 


sgl. 


O9.5-B0.5 lb 


09.7 HI 


n/a 


n/a 


n/a 


CPD -62°2198 


sgl. 


09.5 lb 


09.7 HI 


n/a 


n/a 


n/a 



Notes: a. For comparison purpose, we have inverted the primary and secondary identification so that the mass ratio 
is smaller than unity in the present table. 



to IC 2944 is 2.3 ± 0.3 kpc, in perfect agreement with the 



earlier estimates of Thackeray & Wesselink ( 1965 1 and Tov- 



massian et al. ( 1998 1, but relatively larger than the 1.8 kpc 
estimate of |Kharchenko et al.| ( [2"005[ ). More interestingly, we 
cannot identify significant differences between the various 
subgroups of stars in IC 2944, suggesting that most of the 
O stars belong to a single entity. 

Finally, the dynamical properties of the objects com- 
pared to the average properties of the sample can also be 
used as a membership criterion. From Fig. |19| and beside 
the binaries, only CPD— 62° 2198 displays a significant RV 
difference compared to the systemic velocity of the cluster. 
Combined with the significantly different reddening and dis- 
tance for the object (Table[9|, CPD-62°2198 potentially is 
a much more distant star seen along the line of sight to IC 
2944. 



5.2 Binary fraction 

Eight objects out of 14 O-type stars show a clear SB sig- 
nature, yielding a minimal binary fraction of /min = 0.57. 
Statistical uncertainties due to the limited size of our sam- 
ple amounts to ±0.13 {la error bar). Using the results of 
Sect. |4j the observed binary fraction likely corresponds to a 
true binary fraction /true = 0.65. Proportionally, this makes 
this sample one of the richest binary nest after NGC 6231 
Paper I). 



(Sana et al. 2008 , Paper I). Yet, some objects might 
not belong to IC 2944 itself. For example, HD 100099 and 
HD 101545AB are significantly offset from the cluster core 
and CPD -62°2198 seems much more distant than IC 2944. 
Rejecting these three objects, the SB fraction in IC 2944 
would increase to 0.70. 



Beside HD 101205, HD 101413 and HD 101545, pa^ 



son et al. ( 1998 1 note single star results from speckle obser- 



vations of six targets in Table |8] HD 100099, HD 101131, 
HD 101123, HD 101298, HD 101413, HD 101436. Within a 
radius of 3", the average number of companion per primary 
O star is thus of 0.84 and is again very similar to the value 
obtained in NGC 6231. As in the previous papers of this 
series, close to 70% of the O star population is to be found 
in spectroscopic binary systems. 



5.3 Orbital parameters 

Among the 8 spectroscopic binaries in our sample, three 
have periods up to 10 days, two others have periods of sev- 
eral tens of days, and the last three systems are longer pe- 
riod binaries with, most likely, an orbital time scale of the 
order of months to years. While the sample is likely too 
small to significantly constrain the overall period distribu- 
tion, we note that there is no obvious contradiction with the 
period distribution derived for the binaries in NGC 6231 



(Paper 1 1, neither with the preliminary results of Sana & 
|Evans| ( |2010[ ) . Similarly, the mass-ratio distribution is poorly 
sampled but would qualitatively follow the uniform distri- 



bution suggested by Sana & Evans (20101. 



6 SUMMARY 

Using an extended set of high-resolution high SNR optical 
spectra, we revisited the status of 14 early-type objects in 
IC 2944 and the Cen OB2 association. We presented new 
evidence of binarity for five objects and we confirmed the 
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multiple nature of another two. We derived the first or- 
bital solutions for HD 100099, HD 101436 and HD 101190 
and provide additional support for a higher multiplicity in 
HD 101205. The minimal spectroscopic binary fraction in 
our sample is 0.57 but rises to 0.63 if one only considers 
the objects located in the inner 12' around HD 101205. Our 
sample is limited and does not provide strong constraints 
on the distribution of the orbital parameters of the binary 
systems. Yet, no sizeable disagreement is observed with re- 
sults from earlier papers in this series, nor with the pre- 



liminary results of Sana & Evans (20101 that combine data 



from several clusters. Using an approach similar to the one 
presented in [Paper H[ we estimated the observational biases 
of our campaign and we showed that the binary detection 
rate is close to 90%, leaving thus little room for undetected 
systems. Using newly derived spectroscopic parallaxes, we 
reassessed the distance to IC 2944. We confirm that, as far 
as the O stars are concerned, the cluster is most likely a 
single entity and that HD101545 and HD100099, offset by 
50-60' from the cluster core, are located at a similar distance 
as IC 2944. 
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